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ABSTRACT 

Aims.  We  present  a  unique  multi-epoch  infrared  interferometric  study  of  the  oxygen-rich  Mira  variable  RR  Aql  in  comparison  to 
radiative  transfer  models  of  the  dust  shell.  We  investigate  flux  and  visibility  spectra  at  8-13  //m  with  the  aim  of  better  understanding 
the  pulsation  mechanism  and  its  connection  to  the  dust  condensation  sequence  and  mass-loss  process. 

Methods.  We  obtained  13  epochs  of  mid-infrared  interferometry  with  the  MIDI  instrument  at  the  VLTI  between  April  2004  and  July 
2007,  covering  minimum  to  pre-maximum  pulsation  phases  (0.45-0.85)  within  four  cycles.  The  data  are  modeled  with  a  radiative 
transfer  model  of  the  dust  shell  where  the  central  stellar  intensity  profile  is  described  by  a  series  of  dust-free  dynamic  model  at¬ 
mospheres  based  on  self-excited  pulsation  models.  We  examined  two  dust  species,  silicate  and  ABO3  grains.  We  performed  model 
simulations  using  variations  in  model  phase  and  dust  shell  parameters  to  investigate  the  expected  variability  of  our  mid-infrared  pho¬ 
tometric  and  interferometric  data. 

Results.  The  observed  visibility  spectra  do  not  show  any  indication  of  variations  as  a  function  of  pulsation  phase  and  cycle.  The 
observed  photometry  spectra  may  indicate  intracycle  and  cycle-to-cycle  variations  at  the  level  of  1-2  standard  deviations.  The  pho¬ 
tometric  and  visibility  spectra  of  RR  Aql  can  be  described  well  by  the  radiative  transfer  model  of  the  dust  shell  that  uses  a  dynamic 
model  atmosphere  describing  the  central  source.  The  best-fitting  model  for  our  average  pulsation  phase  of  <lv  =  0.64  ±  0.15  includes 
the  dynamic  model  atmosphere  M2  In  (Tmodei  =  2550  K)  with  a  photospheric  angular  diameter  of  0Phot  =  7.6  ±  0.6  mas,  and  a  silicate 
dust  shell  with  an  optical  depth  of  tv  =  2.8  ±  0.8,  an  inner  radius  of  Rjn  =  4.1  ±  0.7RphOt,  and  a  power-law  index  of  the  density 
distribution  of  p  =  2.6  ±  0.3.  The  addition  of  an  AI2O3  dust  shell  did  not  improve  the  model  fit.  However,  our  model  simulations 
indicate  that  the  presence  of  an  inner  ABO3  dust  shell  with  lower  optical  depth  than  for  the  silicate  dust  shell  can  not  be  excluded. 
The  photospheric  angular  diameter  corresponds  to  a  radius  of  Rphot  =  520/2^  Re  and  an  effective  temperature  of  reff  ~  2420  ±  200  K. 
Our  modeling  simulations  confirm  that  significant  intracycle  and  cycle-to-cycle  visibility  variations  are  not  expected  for  RR  Aql  at 
mid-infrared  wavelengths  within  our  uncertainties. 

Conclusions.  We  conclude  that  our  RR  Aql  data  can  be  described  by  a  pulsating  atmosphere  surrounded  by  a  silicate  dust  shell.  The 
effects  of  the  pulsation  on  the  mid-infrared  flux  and  visibility  values  are  expected  to  be  less  than  about  25%  and  20%,  respectively, 
and  are  too  low  to  be  detected  within  our  measurement  uncertainties. 

Key  words,  techniques:  interferometric  -  stars:  AGB  and  post-AGB  -  stars:  atmospheres  -  stars:  mass-loss  - 
stars:  individual:  RR  Aql 


1.  Introduction 

Asymptotic  giant  branch  (AGB)  stars  are  low-to-intermediate 
mass  stars  at  the  end  of  their  stellar  evolution.  These  stars,  in¬ 
cluding  Mira  type  stars,  exhibit  many  complex  processes  such  as 
shock  fronts  propagating  through  the  stellar  atmosphere,  large- 
amplitude  pulsation,  and  molecule  and  dust  formation  leading  to 
strong  mass  loss  via  a  dense  and  dusty  outflow  from  an  extended 
stellar  atmosphere  (Andersen  et  al.  2003).  The  mass  loss  rates 


*  Based  on  observations  made  with  Very  Large  Telescope 
Interferometer  (VLTI)  at  the  Paranal  Observatory  under  program  IDs 
073.D-071 1,  075.D-0097,  077.D-0630,  and  079.D-0172. 

**  Figure  5  is  available  in  electronic  form  at 
http : / /www . aanda . or g 


can  reach  up  to  1 0  4  MQ  yr  1  (Matsuura  et  al.  2009)  with  expan¬ 
sion  velocities  of  5-30  km  s'1  (Hofner  &  Andersen  2007).  Dust 
grains  are  created  by  condensation  from  the  gas  phase  in  the 
very  cool  and  dense  outer  part  of  the  pulsating  atmosphere.  Dust 
created  by  carbon-rich  stars  owing  to  its  high  opacity  absorbs 
radiation  from  the  star,  and  the  wind  is  believed  to  be  driven  by 
radiation  pressure  on  the  dust  particles  that  drag  gas  along  via 
collisions  (Hofner  2008). 

Despite  remarkable  progress  in  theoretical  and  observational 
studies,  it  still  remains  unclear  whether  the  dust  in  oxygen-rich 
envelopes  is  opaque  enough  to  drive  the  observed  massive  out¬ 
flows  (Woitke  2006;  Hofner  2008).  Because  of  the  wind,  all 
the  matter  around  the  core  of  the  star  is  eventually  returned 
to  the  interstellar  medium  (ISM),  and  the  star  evolves  toward 
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the  planetary  nebula  phase,  leaving  the  core  as  a  white  dwarf. 
Generally  many  aspects  of  the  physics  of  AGB  stars  remain 
poorly  understood,  such  as  the  detailed  atmospheric  stratifica¬ 
tion  and  composition  of  the  stars,  or  the  role  of  stellar  pulsation 
and  its  connection  to  the  dust  formation  and  massive  outflow. 
Improving  our  understanding  of  the  physical  processes  leading 
to  the  substantial  mass  loss  is  important,  as  AGB  stars  play  a  cru¬ 
cial  role  in  the  chemical  enrichment  and  evolution  of  galaxies  by 
returning  gas  and  dust  to  the  ISM. 

Thanks  to  their  large  diameters  and  high  luminosities, 
Mira  variables  are  ideal  targets  for  high  angular  resolu¬ 
tion  observations.  Near-infrared  interferometric  techniques  pro¬ 
vide  detailed  information  regarding  the  conditions  near  the 
continuum-forming  photosphere  such  as  effective  temperature, 
center-to-limb  intensity  variations,  the  stellar  photospheric  di¬ 
ameter,  and  its  dependence  on  wavelength  and  pulsation  phase 
(e.g.,  Haniff  et  al.  1995;  Perrin  et  al.  1999;  van  Belle  et  al. 
1996;  Young  et  al.  2000a;  Thompson  et  al.  2002a, b;  Quirrenbach 
et  al.  1992;  Fedele  et  al.  2005;  Ohnaka  2004;  Millan-Gabet  et  al. 
2005;  Woodruff  et  al.  2008).  New  theoretical  self-excited  dy¬ 
namic  model  atmospheres  of  oxygen-rich  stars  have  been  cre¬ 
ated  and  successfully  applied  (Tej  et  al.  2003;  Hofmann  et  al. 
1998;  Ireland  et  al.  2004b, a,  2008;  Wittkowski  et  al.  2008). 

Observations  at  mid-infrared  wavelengths  are  well-suited  to 
studying  the  molecular  shells  and  the  dust  formation  zone  of 
evolved  stars.  Mid-infrared  interferometry  is  sensitive  to  the 
chemical  composition  and  geometry  of  dust  shells,  their  tem¬ 
perature,  inner  radii,  radial  distribution,  and  the  mass  loss  rate 
(e.g.,  Danchi  et  al.  1994;  Monnier  et  al.  1997, 2000;  Weiner  et  al. 
2006).  Lopez  et  al.  (1997)  presented  long-term  observations  at 
1 1  jim  of  o  Ceti  obtained  with  the  Infrared  Spatial  Interferometer 
(ISI).  The  observed  visibilities  change  from  one  epoch  to  the 
next  and  are  not  consistent  with  simple  heating  or  cooling  of  the 
dust  with  change  in  luminosity  as  a  function  of  stellar  phase, 
but  rather  with  large  temporal  variations  in  the  density  of  the 
dust  shell.  The  data  were  compared  to  axially  symmetric  ra¬ 
diative  transfer  models  and  suggest  inhomogeneities  or  clumps. 
Tevousjan  et  al.  (2004)  has  studied  the  spatial  distribution  of  dust 
around  four  late-type  stars  with  ISI  at  11.15  /rm,  and  find  that  the 
visibility  curves  change  with  the  pulsation  phase  of  the  star.  The 
dust  grains  were  modeled  as  a  mixture  of  silicates  and  graphite. 
The  results  suggest  that  the  dust  shells  appear  to  be  closer  to  the 
star  at  minimum  pulsation  phase,  and  farther  away  at  maximum 
phase,  which  is  also  demonstrated  by  Wittkowski  et  al.  (2007). 
Ohnaka  et  al.  (2007)  reported  on  temporal  visibility  variations  of 
the  carbon-rich  Mira  variable  V  Oph  with  the  instruments  VINCI 
and  MIDI  at  the  VLTI.  This  temporal  variation  of  the  N-band 
angular  sizes  is  largely  governed  by  the  variations  in  the  opacity 
and  the  geometrical  extension  of  the  molecular  layers  (C2H2  and 
HCN)  and  the  dust  shell  (amC  +  SiC). 

Additional  information  to  the  near-  and  mid-infrared  obser¬ 
vations  can  be  obtained  with  Very  Long  Baseline  Array  (VLBA) 
observations,  which  allow  one  to  study  the  properties  of  the  cir- 
cumstellar  environment  of  evolved  stars  using  the  maser  radia¬ 
tion  emitted  by  some  molecules,  most  commonly  SiO,  H2O,  and 
OH.  By  observing  the  maser  radiation  at  tens  to  several  hundred 
AU  from  the  star,  it  is  possible  to  determine  accurate  proper  mo¬ 
tions  and  parallaxes  (Vlemmings  et  al.  2003),  total  intensity  and 
linear  polarizations  (Cotton  et  al.  2008),  and  the  structure  of  the 
environment  of  the  stars  (Diamond  et  al.  1994). 

This  work  is  part  of  the  ongoing  project  of  concurrent  mul¬ 
tiwavelength  observations  using  optical  long  baseline  interfer¬ 
ometry  (AMBER  and  MIDI  at  the  VLTI)  combined  with  ra¬ 
dio  interferometry  (VLBA).  We  aim  to  investigate  layers  at 
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Fig.  l.  Visual  light  curve  of  RR  Aql  based  on  data  from  the  AAVSO  and 
AFOEV  databases  as  a  function  of  Julian  Date  and  stellar  cycle/phase. 
The  green  arrows  indicate  the  dates  of  our  VLTI/MIDI  observations. 
The  blue  arrows  indicate  the  dates  of  our  VLBA  observations.  Here, 
the  full  arrows  denote  observations  of  SiO  maser  emission  and  simple 
arrows  observations  of  H2O  maser  emission. 


different  depths  of  the  atmosphere  and  circumstellar  envelope 
(CSE)  of  the  stars.  Previous  results  from  this  project  of  joint 
VLTI/VLBA  observations  for  the  Mira  star  S  Ori  can  be  found 
in  Boboltz  &  Wittkowski  (2005)  and  Wittkowski  et  al.  (2007). 
Here  we  present  results  of  a  long-term  VLTI/MIDI  monitor¬ 
ing  of  the  oxygen-rich  Mira  variable  RR  Aql.  In  addition  to 
the  VLTI/MIDI  observation  presented  in  this  paper,  we  obtained 
coordinated  multi-epoch  VLBA  observations  of  SiO  and  H2O 
masers  towards  RR  Aql,  which  will  be  presented  in  a  subsequent 
paper. 


2.  Characteristics  of  RR  Aql 

RR  Aql  is  an  oxygen-rich  Mira  variable  with  spectral  type  M6e- 
M9  (Samus  et  al.  2004).  RR  Aql  shows  a  strong  silicate  emis¬ 
sion  feature  in  its  mid-infrared  spectrum  (Lorenz-Martins  & 
Pompeia  2000).  In  addition,  it  has  relatively  strong  SiO,  H2O, 
and  OH  maser  emission  (Benson  et  al.  1990).  Vlemmings  &  van 
Langevelde  (2007)  observed  OH  masers  toward  RR  Aql  at  five 
epochs  with  the  VLBA.  Based  on  these  observations,  the  dis¬ 
tance  to  RR  Aql  is  estimated  to  D  -  633  pc  (HIPPARCOS 
distance  D  -  540  pc).  Van  Belle  et  al.  (2002)  obtained  a  JGband 
(A  =  2.2jitm,  AA  —  0.4 /mi)  angular  size  with  the  Infrared 
Optical  Telescope  Array  (IOTA).  The  uniform  disk  diameter 
(UD)  of  0ud  =  10.73  ±  0.66  mas  at  phase  ®  =  0.48,  an  ef¬ 
fective  temperature  Teg  =  2127  ±  111K,  and  a  bolometric  flux 
/bol  =  78.4  ±  11.8  x  10~8ergs  cm'V.  Miyata  et  al.  (2000) 
studied  the  dust  around  the  star  spectroscopically  and  obtained 
Edust/Fstar  at  10/rm  of  1.49  ±  0.02.  Ragland  et  al.  (2006)  mea¬ 
sured  the  closure  phase  with  IOTA  in  the  H- band,  and  classified 
RR  Aql  as  a  target  with  no  detectable  asymmetries.  The  IRAS 
flux  at  1 2  fjm  is  332  Jy.  The  light  curve  in  the  V  band  varies 
from  ~8  mag  at  maximum  light  to  ~  14  mag  at  minimum  light. 
Whitelock  et  al.  (2000)  shows  that  the  K  magnitude  of  RR  Aql 
varies  from  ~0.0  mag  at  maximum  light  to  ~1.0  mag  at  mini¬ 
mum  light.  RR  Aql  is  pulsating  with  a  period  of  P  =  394.78  days 
(Samus  et  al.  2004)  and  the  Julian  Date  of  the  last  maximum 
brightness  is  Tq  =  2452875.4  (Pojmanski  et  al.  2005).  The  pe¬ 
riod  of  pulsation  visually  corresponds  to  recent  values  from  the 


A134,  page  2  of  15 


I.  Karovicova  et  al.:  Mid-infrared  interferometric  monitoring  of  the  Mira  variable  RR  Aql 


Table  1.  VLTI/MIDI  observation  of  RR  Aql. 


Ep. 

date 

DD/MM/YYYY 

Time 

[UTC] 

JD 

®vis 

Conf. 

B 

Disp. 

Elem. 

BC 

BP 

[m] 

PA 

[deg] 

Seeing 

["] 

To 

[ms] 

A 

09/04/2004 

07:26 

2453105 

0.58 

U2-U3 

Aim 

Prism 

HS 

33.57 

3.21 

0.45 

7.4 

A 

10/04/2004 

09:05 

2453106 

0.58 

U2-U3 

47m 

Prism 

HS 

37.17 

26.09 

0.76 

7.8 

A 

09:52 

2453106 

0.58 

U2-U3 

47m 

Prism 

HS 

39.86 

33.48 

0.64 

8.2 

B 

09/07/2004 

08:25 

2453196 

0.81 

U2-U3 

47m 

Prism 

HS 

45.23 

44.83 

0.47 

4.8 

C 

28/07/2004 

08:09 

2453215 

0.86 

U2-U3 

47m 

Prism 

HS 

42.54 

41.57 

1.36 

1.0 

C 

29/07/2004 

02:05 

2453216 

0.86 

U2-U3 

47m 

Prism 

HS 

37.83 

28.16 

1.66 

1.0 

C 

03:33 

2453216 

0.86 

U2-U3 

47m 

Prism 

HS 

42.84 

39.49 

0.75 

2.3 

C 

04:55 

2453216 

0.86 

U2-U3 

47m 

Prism 

HS 

45.99 

44.67 

1.04 

1.6 

C 

06:33 

2453216 

0.86 

U2-U3 

47m 

Prism 

HS 

46.16 

45.71 

3.38 

0.5 

C 

07:04 

2453216 

0.86 

U2-U3 

47m 

Prism 

HS 

45.32 

44.92 

2.68 

0.7 

C 

01/08/2004 

00:12 

2453219 

0.87 

U2-U3 

47m 

Prism 

HS 

33.76 

6.86 

0.73 

3.1 

D 

18/04/2005 

10:30 

2453479 

1.53 

U2-U4 

89m 

Prism 

HS 

88.37 

81.72 

0.76 

3.5 

D 

19/04/2005 

07:53 

2453480 

1.53 

U2-U4 

89m 

Prism 

HS 

61.36 

76.41 

0.83 

2.3 

E 

20/07/2005 

03:51 

2453572 

1.76 

U1-U4 

130m 

Prism 

HS 

119.63 

59.78 

0.71 

1.9 

E 

06:16 

2453572 

1.76 

U1-U4 

130m 

Prism 

HS 

128.64 

63.28 

0.81 

1.6 

E 

22/07/2005 

05:10 

2453574 

1.77 

U2-U3 

47m 

Prism 

HS 

45.65 

44.12 

1.00 

1.4 

E 

07:56 

2453574 

1.77 

U2-U3 

47m 

Prism 

HS 

44.44 

43.96 

1.11 

1.2 

F 

18/04/2006 

09:38 

2453844 

2.45 

D0-G0 

32m 

Prism 

HS 

28.99 

70.01 

1.08 

2.1 

F 

19/04/2006 

08:32 

2453845 

2.46 

D0-G0 

32m 

Prism 

HS 

24.74 

65.81 

0.80 

6.0 

F 

09:21 

2453845 

2.46 

D0-G0 

32m 

Prism 

HS 

28.21 

69.32 

0.78 

6.3 

G 

21/05/2006 

05:53 

2453877 

2.54 

E0-G0 

16m 

Prism 

HS 

10.94 

62.10 

0.56 

4.0 

G 

25/05/2006 

06:35 

2453881 

2.55 

A0-G0 

64m 

Prism 

HS 

53.17 

67.77 

0.55 

5.1 

G 

08:46 

2453881 

2.55 

A0-G0 

64m 

Prism 

HS 

63.93 

72.67 

0.44 

6.7 

H 

18/06/2006 

03:58 

2453905 

2.61 

D0-G0 

32m 

Prism 

HS 

21.37 

61.34 

1.29 

1.5 

H 

04:41 

2453905 

2.61 

D0-G0 

32m 

Prism 

HS 

25.09 

66.21 

1.27 

1.6 

H 

18/06/2006 

05:39 

2453905 

2.61 

D0-G0 

32m 

Prism 

HS 

29.05 

70.06 

1.00 

2.1 

H 

20/06/2006 

05:49 

2453907 

2.61 

A0-G0 

64m 

Prism 

HS 

59.89 

70.81 

1.05 

3.9 

H 

06:33 

2453907 

2.61 

A0-G0 

64m 

Prism 

HS 

63.00 

72.15 

0.88 

4.5 

H 

21/06/2006 

07:08 

2453908 

2.62 

E0-G0 

16m 

Prism 

HS 

16.00 

72.76 

1.01 

2.2 

H 

09:36 

2453908 

2.62 

E0-G0 

16m 

Prism 

HS 

13.21 

70.37 

0.77 

2.9 

H 

23/06/2006 

03:20 

2453910 

2.62 

E0-G0 

16m 

Prism 

HS 

9.84 

58.42 

0.76 

3.9 

H 

04:07 

2453910 

2.62 

E0-G0 

16m 

Prism 

HS 

11.97 

64.86 

0.67 

4.4 

I 

08/08/2006 

05:12 

2453956 

2.74 

A0-G0 

64m 

Prism 

HS 

61.41 

72.62 

1.02 

2.3 

I 

05:57 

2453956 

2.74 

A0-G0 

64m 

Prism 

HS 

56.88 

71.58 

1.10 

2.2 

I 

09/08/2006 

06:39 

2453957 

2.74 

A0-G0 

64m 

Prism 

HS 

50.33 

69.49 

1.47 

1.5 

I 

10/08/2006 

03:07 

2453958 

2.74 

DO-GO 

32m 

Prism 

HS 

31.37 

72.03 

1.43 

2.6 

I 

04:31 

2453958 

2.74 

DO-GO 

32m 

Prism 

HS 

31.65 

72.90 

1.67 

2.5 

I 

05:45 

2453958 

2.74 

DO-GO 

32m 

Prism 

HS 

28.72 

71.74 

1.43 

3.2 

I 

1 1/08/2006 

02:44 

2453959 

2.74 

DO-GO 

32m 

Prism 

HS 

30.76 

71.51 

1.43 

2.3 

I 

13/08/2006 

04:05 

2453961 

2.75 

EO-GO 

16m 

Prism 

HS 

15.94 

72.90 

0.91 

1.9 

I 

04:56 

2453961 

2.75 

EO-GO 

16m 

Prism 

HS 

15.26 

72.56 

0.63 

2.7 

I 

05:52 

2453961 

2.75 

EO-GO 

16m 

Prism 

HS 

13.76 

71.07 

0.60 

2.8 

J 

14/09/2006 

01:01 

2453993 

2.83 

EO-GO 

16m 

Prism 

HS 

15.82 

72.29 

0.85 

2.3 

J 

16/09/2006 

02:15 

2453995 

2.84 

AO-GO 

64m 

Prism 

HS 

62.86 

72.85 

0.95 

1.6 

J 

03:13 

2453995 

2.84 

AO-GO 

64m 

Prism 

HS 

58.12 

71.90 

1.68 

0.9 

K 

09/05/2007 

08:53 

2454230 

3.43 

DO-HO 

64m 

Prism 

HS 

61.43 

71.46 

1.12 

1.3 

L 

22/06/2007 

03:50 

2454274 

3.54 

GO-HO 

32m 

Prism 

HS 

21.96 

62.20 

1.17 

1.0 

L 

05:29 

2454274 

3.54 

GO-HO 

32m 

Prism 

HS 

29.30 

70.26 

1.21 

0.9 

M 

03/07/2007 

07:25 

2454285 

3.57 

EO-GO 

16m 

Prism 

HS 

15.51 

72.73 

1.07 

2.8 

M 

07:37 

2454285 

3.57 

EO-GO 

16m 

Prism 

HS 

15.31 

72.59 

0.75 

3.9 

M 

04/07/2007 

04:42 

2454286 

3.57 

GO-HO 

32m 

Prism 

HS 

29.28 

70.24 

1.59 

1.0 

M 

06:30 

2454286 

3.57 

GO-HO 

32m 

Prism 

HS 

31.98 

72.83 

1.51 

1.1 

Notes.  The  table  lists  the  epoch,  the  date,  the  time,  the  Julian  Date  (JD),  the  visual  pulsation  phase  ®vis,  the  baseline  configuration,  the  ground 
length  of  the  configuration,  the  dispersive  element,  the  beam  combiner  BC,  the  projected  baseline  length  6p,  the  position  angle  on  the  sky  PA  (deg. 
east  of  north),  the  DIMM  seeing  (at  500  nm),  and  the  coherence  time  r0  (at  500  nm). 


AAVSO1  and  AFOEV2  databases.  Figure  1  shows  the  visual 
light  curve  of  RR  Aql  based  on  values  from  the  AAVSO  and 
AFOEV  databases  as  a  function  of  Julian  day  and  stellar  phase. 


1  http://www.aavso.org 

2  http : //cdsweb . u- strasbg . fr/ afoev 


3.  VLTI/MIDI  observations  and  data  reduction 

We  obtained  57  spectrally-dispersed  mid-infrared  interferomet¬ 
ric  observations  of  RR  Aql  with  the  VLTI/MIDI  instrument  be¬ 
tween  Apr.  9,  2004  and  Jul.  28,  2007.  MIDI  -  the  mid-infrared 
interferometric  instrument  (Leinert  et  al.  2004)  combines  the 
beams  from  two  telescopes  of  the  VLTI  (Glindemann  et  al.  2003) 
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Fig.  2.  Calibrated  RR  Aql  MIDI  flux  spectra  as  a  function  of  wave¬ 
length.  For  clarity,  the  error  bars  are  omitted  in  the  plot. 


Wavelength  [jU.ni] 


Wavelength  [jum] 

Fig.  3.  RR  Aql  flux  spectrum  as  a  function  of  wavelength  from  2.4  p m  to 
40yum  {top).  The  lines  represent  the  flux  spectra  from  ISO-SWS  obser¬ 
vations  covering  wavelengths  from  2.5  pm  to  ~240 pm  (solid  red  thick 
line),  IRAS-LRS  observations  covering  ~7 .7 pm  to  ~23pm  (dashed 
thin  black  line),  and  the  mean  of  our  (V-band  MIDI  measurements  (solid 
thin  green  line).  Here,  the  vertical  bars  span  the  maximum  and  mini¬ 
mum  values  measured.  The  diamonds  denote  2MASS  measurements  at 
1.25  pm,  1.65  pm,  and  2.2  pm.  The  solid  black  line  indicates  our  atmo¬ 
sphere  and  dust  shell  model  as  explained  in  Sect.  5.2.  The  bottom  plot 
shows  an  enlarged  segment  of  the  plot  in  the  MIDI  wavelength  range  of 
8-13  pm. 


and  provides  spectrally  resolved  visibilities  in  the  10  pm  window 
(/V-band,  8-1 3 pm).  To  obtain  dispersed  photometric  and  inter¬ 
ferometric  signals,  we  used  the  PRISM  as  a  dispersive  element 
with  a  spectral  resolution  R  —  A  A/ A  ~  30.  The  beams  were 
combined  in  High_Sens  mode  (HS).  In  this  mode  the  photomet¬ 
ric  signal  is  observed  after  the  interferometric  signal. 


The  details  of  the  observations  and  the  instrumental  settings 
are  summarized  in  Table  1.  All  observations  were  executed  in 
service  mode  using  either  the  unit  telescopes  (UTs,  8.2  m)  or  the 
auxiliary  telescopes  (ATs,  1.8  m). 

We  merged  the  MIDI  data  into  15  epochs,  with  a  maximum 
time-lag  between  individual  observations  of  five  days  for  each 
epoch  ( 1 .3%  of  the  pulsation  period).  For  technical  problems  we 
had  to  eliminate  two  epochs.  Figure  1  shows  the  final  13  epochs 
in  comparison  to  the  light  curve.  The  uncertainty  in  the  allo¬ 
cation  of  the  visual  phase  to  our  observations  was  estimated  to 
~0.1.  Our  long-term  VLBA  monitoring  of  SiO  and  H2O  maser 
emission  toward  RR  Aql,  will  be  presented  in  a  forthcoming 
paper. 

We  used  the  MIA+EWS  software  package,  version  1.6  (Jaffe 
et  al.3)  for  the  MIDI  data  reduction.  This  package  includes 
two  different  methods,  an  incoherent  method  (MPIA  software 
package  MIA)  that  analyzes  the  powerspectrum  of  the  observed 
fringe  signal  and  a  coherent  integration  method  (EWS),  which 
first  compensates  for  optical  path  differences,  including  both  in¬ 
strumental  and  atmospheric  delays  in  each  scan,  and  then  co¬ 
herently  adds  the  fringes.  We  applied  both  methods  to  indepen¬ 
dently  verify  the  data  reduction  results.  The  detector  masks  were 
calculated  by  the  procedure  of  MIA,  and  were  used  for  both  the 
MIA  and  the  EWS  analysis.  The  obtained  data  reduction  values 
correspond  to  each  other,  and  we  chose  to  use  in  the  follow¬ 
ing  the  results  derived  from  the  EWS  analysis,  which  offer  error 
estimations. 

To  account  for  instrumental  visibility  losses  and  to  deter¬ 
mine  the  absolute  flux  values,  calibrator  stars  with  known  flux 
and  diameter  values  were  observed  immediately  before  or  af¬ 
ter  the  science  target.  Our  main  calibrators  were  HD  169916 
(period  P73,  P75),  HD  146051  (P77),  and  HD  177716  (P79). 
Calibrated  science  target  visibility  spectra  were  calculated  us¬ 
ing  the  instrumental  transfer  function  derived  from  all  calibrator 
data  sets  taken  during  the  same  night  with  the  same  baseline  and 
instrumental  mode  as  our  scientific  target.  The  number  of  avail¬ 
able  transfer  function  measurements  depends  on  the  number  of 
calibrator  stars  observed  per  specific  night  including  those  cal¬ 
ibrators  observed  by  other  programs.  The  errors  of  the  transfer 
functions  are  given  by  the  standard  deviation  of  all  transfer  func¬ 
tion  measurements  per  night.  To  estimate  the  uncertainty  of  the 
transfer  function  for  nights  when  only  one  calibrator  was  avail¬ 
able,  we  used  typical  values  based  on  nights  when  many  calibra¬ 
tor  stars  were  observed.  The  final  errors  on  the  observed  visibili¬ 
ties  are  mostly  systematic,  and  include  the  error  of  the  coherence 
factor  of  the  science  target  and  the  calibrators,  the  adopted  di¬ 
ameter  errors,  and  the  standard  deviation  of  the  transfer  function 
over  the  night. 

The  photometric  spectrum  was  calibrated  with  one  or  two 
calibrators,  which  were  observed  close  on  the  sky  and  in  time 
compared  to  our  science  target.  For  most  calibrator  stars,  abso¬ 
lutely  calibrated  spectra  are  available  in  Cohen  et  al.  (1999).  For 
those  calibrators  where  the  absolutely  calibrated  spectrum  was 
not  directly  available,  we  instead  used  a  spectrum  of  a  calibra¬ 
tor  with  a  similar  spectral  type  and  similar  effective  temperature 
(see  the  instrument  consortium’s  catalog4).  The  spectra  of  such 
calibrators  were  scaled  with  the  IRAS  flux  at  1 2  pm  to  the  level 
of  our  calibrator.  In  addition,  we  verified  that  our  synthetic  spec¬ 
tra  obtained  by  this  procedure  are  valid  by  scaling  known  spectra 
of  two  Cohen  calibrators.  In  a  few  cases,  when  the  atmospheric 


3  http: //www. strw. leidenuniv.nl/~nevec/MIDI 

4  http : / /www . ster . kuleuven . ac . be/~ti j 1/MIDI_ 
calibration/mcc . txt 
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Fig.  4.  Calibrated  RR  Aql  MIDI  visibility  amplitudes  as  a  function  of  wavelength.  ( top  left)  Observations  executed  at  projected  baselines  (Bp)  < 
20  m,  ( top  right)  20  m  <  Bp  <  35  m,  ( bottom  left)  35  m  <  Bp  <  50  m,  and  ( bottom  right)  Bp  >  50  m. 


absorption  was  strongly  affecting  the  spectra  around  9.5  /im,  we 
used  another  similar  calibrator  instead  of  the  main  calibrator  ob¬ 
served  in  the  same  night  with  the  same  level  of  flux.  The  ambi¬ 
ent  conditions  for  all  the  observations  were  carefully  checked.  If 
we  were  facing  problems  with  clouds,  constraints  due  to  wind, 
significant  differences  in  seeing,  humidity,  coherence  time,  and 
airmass  between  the  science  target  and  the  corresponding  cali¬ 
brator,  the  photometry  was  omitted  from  the  analysis. 

4.  MIDI  results 

We  observed  RR  Aql  at  different  phases  over  four  pulsation  cy¬ 
cles  in  order  to  monitor  the  mid-infrared  photometry  and  visibil¬ 
ity  spectra.  Here,  we  present  and  discuss  the  general  properties 
of  the  data,  followed  by  an  analysis  of  their  variability  as  a  func¬ 
tion  of  phase  and  cycle  in  Sects.  4.1  and  4.2. 

Figure  2  shows  all  obtained  calibrated  photometry  spectra  as 
a  function  of  wavelength.  The  MIDI  flux  measurements  show  a 
consistent  shape  exhibiting  an  increase  in  the  flux  from  ~  100- 
200  Jy  at  8  /mi  to  a  maximum  near  9.8  pm  of  ~200-400  Jy,  and  a 
decrease  towards  13  /mi,  where  the  flux  values  again  reach  val¬ 
ues  of  ~  1 00-200  Jy.  The  level  of  the  flux  spectrum  differs  for  in¬ 
dividual  measurements  with  a  spread  of  ~1 00-200  Jy.  Figure  3 
shows  a  comparison  of  the  mean  of  our  MIDI  flux  measure¬ 
ments  to  measurements  obtained  with  the  ISO  and  IRAS  in¬ 
struments.  The  shape  of  the  flux  curve  is  consistent  among  the 
MIDI,  ISO,  and  IRAS  measurements.  The  level  of  the  IRAS  flux 
is  within  the  range  of  our  MIDI  measurements.  The  level  of  the 
ISO  flux  is  higher,  which  can  most  likely  be  explained  by  the 
post-maximum  phase  of  0. 16  of  the  ISO  observations  (1997-05- 
03)  compared  to  our  minimum  to  pre-maximum  phases.  The  dif¬ 
ference  between  the  MIDI  and  ISO  flux  level  may  also  indicate 


a  flux  variation  over  different  cycles.  Figure  3  also  includes  a 
model  description  of  our  MIDI  data,  which  is  explained  below 
in  Sect.  5.2. 

Figure  4  shows  all  obtained  calibrated  visibility  spectra, 
which  are  combined  into  four  groups  of  different  projected  base¬ 
line  lengths  (Bp)  of  Bp  <  20  m,  20  m  <  Bp  <  35  m,  35  m  < 
Bp  <  50  m,  and  Bp  >  50  m.  The  visibility  curves  show  a  signifi¬ 
cant  wavelength  dependence  with  a  steep  decrease  from  8  pm  to 
~9.5  pi n  and  a  slow  increase  in  the  9.5  pm  to  13  pm  range.  The 
shape  and  absolute  scale  of  the  visibility  function  depends  on  the 
length  of  the  projected  baseline  where  a  longer  baseline  results 
in  a  lower  and  flatter  curve. 

As  a  first  basic  interpretation  of  the  interferometric  data,  we 
computed  the  corresponding  uniform  disk  (UD)  diameter  and  the 
Gaussian  FWHM  for  each  data  set  and  spectral  channel.  This 
leads  to  a  rough  estimate  of  the  characteristic  size  of  the  tar¬ 
get  at  each  wavelength.  However,  it  should  be  mentioned  that 
the  true  intensity  distribution  across  the  stellar  disk  is  expected 
to  be  more  complex  than  can  be  described  by  these  elemen¬ 
tary  models.  Figure  5  shows  for  the  example  of  epoch  A,  and 
epoch  I  the  flux,  the  visibility  amplitude,  the  corresponding  UD 
diameter,  and  the  corresponding  Gaussian  FWHM  diameter  as 
a  function  of  wavelength.  The  shape  of  the  UD  diameter  and 
Gaussian  FWHM  functions  show  a  steep  increase  from  8  pm  to 
~9.5  pm  by  a  factor  of  2.  A  plateau  appears  between  ~9.5  pm 
and  11.5  pm,  and  the  values  are  nearly  constant  from  ~  1 .5 pm  to 
13  pm. 

The  shape  of  the  8- 13  pm  flux  spectrum  with  a  maximum 
near  9.8pm  is  known  to  be  a  characteristic  silicate  emission 
feature  (e.g.,  Little-Marenin  &  Little  1990;  Lorenz-Martins  & 
Pompeia  2000).  The  corresponding  drop  in  the  visibility  func¬ 
tion  between  8  pm  and  9.5  pm  and  increase  in  the  UD  diameter 
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Fig.  6.  Calibrated  MIDI  visibility  amplitudes  for  different  pulsation 
phases  within  the  same  cycle  to  investigate  intracycle  visibility  varia¬ 
tions.  Each  line  represents  a  different  pulsation  phase  within  the  same 
cycle  and  is  computed  as  an  average  of  data  obtained  at  the  respective 
phase  (±0.15)  and  observed  at  similar  projected  baseline  length  ( Bv  ± 
10%)  and  position  angle  (PA  ±  10%).  The  top  panel  shows  the  exam¬ 
ple  of  pulsation  phases  0.45,  0.61,  and  0.74  of  cycle  two  observed  with 
a  projected  baseline  length  of  ~29m  and  a  position  angle  of  ~70  deg. 
The  bottom  panel  shows  the  example  of  pulsation  phases  0.62,  0.75,  and 
0.83  of  cycle  two  observed  with  a  projected  baseline  length  of  ~16m 
and  a  position  angle  of  ~72  deg.  The  error  bars  are  computed  as  the 
standard  deviation  of  the  averaged  visibilities. 


and  Gaussian  FWHM  has  also  been  shown  to  be  a  typical  sig¬ 
nature  of  a  circumstellar  dust  shell  that  is  dominated  by  silicate 
dust  (cf.,  e.g.,  Driebe  et  al.  2008;  Ohnaka  et  al.  2008).  The  op¬ 
posite  trend  toward  a  broad  spectrum  and  visibility  function  in 
the  wavelength  range  8-13 pm  was  interpreted  as  a  dust  shell 
consisting  of  AI2O3  dust  (Wittkowski  et  al.  2007)  or  silicate  and 
AFO3  dust  (Ohnaka  et  al.  2005). 

4. 1 .  Visibility  monitoring 

We  obtained  a  rich  sample  of  MIDI  data  on  RR  Aql,  which 
covers  a  total  of  four  pulsation  cycles  and  pulsation  phases  be¬ 
tween  0.45  and  0.85,  i.e.,  minimum  to  pre-maximum  phases  (see 
Table  1  and  Fig.  1).  For  many  different  pulsation  cycles  and 
phases,  we  obtained  data  at  similar  projected  baseline  lengths 
and  position  angles.  This  gives  us  a  unique  opportunity  to  mean¬ 
ingfully  compare  interferometric  data  obtained  at  different  pul¬ 
sation  phases  and  cycles.  Since  the  visibility  depends  on  the 
probed  point  on  the  uv  plane,  and  thus  on  the  projected  baseline 
length  and  position  angle,  visibility  data  at  different  pulsation 
phases  can  be  directly  compared  only  if  they  were  obtained  at 
the  same,  or  very  similar,  point  on  the  uv  plane.  For  this  reason. 


Fig.  7.  As  for  Fig.  6,  but  for  the  same  pulsation  phase  in  consecutive 
pulsation  cycles  to  investigate  cycle-to-cycle  visibility  variations.  The 
top  panel  shows  the  example  of  pulsation  phase  ~0.5  in  three  consec¬ 
utive  cycles  observed  with  a  projected  baseline  length  of  ~61  m  and  a 
position  angle  ~75  deg.  The  bottom  panel  shows  the  example  of  pul¬ 
sation  phase  0.5  in  two  consecutive  cycles  observed  with  a  projected 
baseline  length  of  ~29  m  and  a  position  angle  ~70  deg. 


we  combined  individual  observations  into  groups  of  similar  pul¬ 
sation  phases  (®vis+0.15),  projected  baseline  lengths  (Bp±10%), 
and  position  angles  (PA  ±  10%).  The  data  within  each  group 
were  averaged.  The  uncertainty  of  the  averaged  visibility  curves 
was  estimated  as  the  standard  deviation  of  the  averaged  visibili¬ 
ties. 

4.1 .1 .  Intracycle  visibility  monitoring 

To  investigate  intracycle  visibility  variations  we  compared  data 
observed  at  different  pulsation  phases  within  the  same  cycle. 
Figure  6  shows  two  examples  of  calibrated  visibility  curves, 
where  each  line  in  the  plot  represents  an  average  of  a  group  of 
visibility  data,  as  described  above,  (i.e.  with  similar  Bp,  PA.  and 
®vis).  The  top  panel  shows  an  example  of  observations  within 
the  second  pulsation  cycle  at  phases  2.45,  2.61,  and  2.74  ob¬ 
tained  with  Bp  ~  30  m  and  PA  ~  70°.  The  bottom  panel  shows 
an  example  of  observations  at  phases  2.62,  2.75,  and  2.83  ob¬ 
tained  with  Bp  ~  16  m  and  PA  ~  72°.  The  data  do  not  show 
any  evidence  of  intracycle  visibility  variations  within  the  probed 
range  of  pulsation  phases  (~0. 45-0. 85)  and  within  our  visibility 
accuracies  of  about  5-20%. 

4.1.2.  Cycle-to-cycle  visibility  monitoring 

We  compared  data  observed  at  similar  pulsation  phases  of  dif¬ 
ferent  consecutive  cycles  to  investigate  cycle-to-cycle  visibility 
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Fig.  8.  Calibrated  MIDI  flux  spectrum  for  different  pulsation  phases 
within  the  same  cycle  to  investigate  intracycle  photometry  variations. 
Each  line  represents  a  different  pulsation  phase  within  the  same  cycle 
and  is  computed  as  an  average  of  data  obtained  at  the  respective  phase 
(±0.15).  The  top  panel  shows  the  example  of  pulsation  phases  0.58  and 
0.85  of  cycle  zero.  The  bottom  panel  shows  the  examples  of  phases  0.53 
and  0.76  of  cycle  one.  The  error  bars  are  computed  as  the  standard  devi¬ 
ation  of  the  averaged  photometry  curves.  The  gray  shades  denote  zones 
that  are  affected  by  atmospheric  absorption. 


variations.  Figure  7  shows  two  of  these  examples  computed  in 
the  same  way  as  shown  in  Fig.  6  of  Sect.  4.1.1,  but  (top  panel)  for 
minimum  phases  of  1.53,  2.59,  and  3.43  in  three  consecutive  cy¬ 
cles  obtained  with  a  projected  baseline  length  of  ~  62  m  and  PA 
~  73°,  and  (bottom  panel)  for  minimum  phases  2.43,  and  3.55  of 
two  consecutive  cycles  obtained  with  a  projected  baseline  length 
of  -29  m,  and  PA  ~  70°.  Here,  the  phases  of  the  different  cycles 
differ  by  up  to  ~20%.  However,  we  showed  in  Sect.  4.1.1  that 
there  is  no  evidence  of  intracycle  visibility  variations,  so  that  the 
chosen  phases  can  be  compared  well  for  different  cycles.  As  a 
result,  Fig.  7  shows  that  our  data  do  not  exhibit  any  significant 
cycle-to-cycle  visibility  variation  for  minimum  phases  over  two 
or  three  cycles  within  our  visibility  accuracies  of  ~5-20%. 

4.1 .3.  Deviations  from  circular  symmetry 

Most  position  angles  varied  around  40°  or  70°.  Only  in  two 
cases  is  the  PA  significantly  different  with  values  of  3°  and  7°. 
However,  in  both  cases  the  visibility  data  quality  is  poor.  The 
differential  phases  are  close  to  zero  within  ~  10-20  deg,  which 
is  close  to  the  calibration  uncertainties  that  can  be  reached  for 
MIDI  differential  phases  (e.g.,  Ohnaka  et  al.  2008).  Therefore 
these  observations  do  not  allow  us  to  make  any  reliable  conclu¬ 
sion  about  the  presence  or  absence  of  an  asymmetric  intensity 
distribution. 


Wavelength  [/zm] 

Fig.  9.  As  for  Fig.  8,  but  for  the  same  pulsation  phase  in  consecutive  pul¬ 
sation  cycles  to  investigate  cycle-to-cycle  photometry  variations.  The 
top  panel  shows  the  example  of  phase  ~0.75  in  cycles  one  and  two.  The 
bottom  panel  shows  the  example  of  phase  ~0.85  of  cycles  zero  and  two. 


Summarizing,  our  RR  Aql  MIDI  data  do  not  show  any  evi¬ 
dence  of  significant  intracycle  or  cycle-to-cycle  visibility  varia¬ 
tions  within  the  examined  phase  coverage  of  ~0. 45-0. 85  in  four 
consecutive  cycles  (see  Table  1  for  details  of  the  phase  coverage) 
and  within  our  visibility  accuracies  of  ~5-20%.  The  obtained 
mid-infrared  interferometric  observations  imply  either  that  the 
mid-infrared  sizes  of  molecular  and  dusty  layers  of  RR  Aql 
do  not  significantly  vary  within  our  phase  coverage  or  that  the 
conducted  observations  are  not  sensitive  enough  to  detect  such 
variations.  In  addition,  the  good  agreement  of  visibility  data  ob¬ 
tained  at  the  same  projected  baseline  lengths  and  position  angles 
for  several  years  also  confirms  the  good  data  quality  and  credi¬ 
bility  of  the  data  reduction  procedure. 

4.2.  N-band  flux  monitoring 

To  investigate  the  intracycle  and  cycle-to-cycle  variability  of 
our  8-1 3  pm  flux  spectra,  we  combined  the  data  into  groups  of 
data  obtained  at  similar  pulsation  phases  (<Dvjs  ±0.15).  The  data 
within  each  group  were  averaged.  The  uncertainty  of  the  aver¬ 
aged  photometry  curves  was  estimated  as  the  standard  deviation 
of  the  averaged  values. 

4.2.1.  Intracycle  photometry  monitoring 

Figure  8  shows  two  examples  of  a  comparison  of  the  MIDI  cal¬ 
ibrated  flux  spectrum  at  different  pulsation  phases  of  the  same 
cycle.  The  top  panel  shows  a  comparison  between  phases  0.58 
and  0.85  of  cycle  zero,  and  the  bottom  panel  shows  a  comparison 
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between  phases  1.53  and  1.76  of  cycle  one.  Both  examples 
compare  data  taken  at  post-maximum  and  pre-minimum  phases, 
i.e.  with  the  largest  separation  in  phase  that  is  available  in  our 
dataset.  The  RR  Aql  flux  values  within  the  yV-band  are  larger  at 
the  pre-maximum  phases  than  at  the  post-minimum  phases,  cor¬ 
responding  to  the  light  curve  of  the  V-magnitude.  The  difference 
is  most  pronounced  toward  the  silicate  emission  feature  around 
9.8/,/m  and  smaller  towards  the  edges  of  the  MIDI  bandpass  at 
8.0/rm  and  1 3.0  pm.  This  can  most  likely  be  understood  as  a 
consequence  of  the  silicate  absorption  coefficient  being  largest 
at  short  visual  wavelengths,  so  that  the  silicate  temperature  is 
particularly  sensitive  to  the  large  visual  flux  variations  of  the 
central  star.  The  differences  correspond  to  up  to  20%-35%,  or 
about  one  or  two  standard  deviations. 


4.2.2.  Cycle-to-cycle  photometry  monitoring 

Figure  9  shows  a  comparison  of  data  observed  at  the  same  phase 
of  consecutive  pulsation  cycles.  The  top  panel  shows  a  compar¬ 
ison  between  phases  1.76  and  2.74,  i.e.  phase  0.75  of  cycles  one 
and  two.  The  bottom  panel  shows  a  comparison  between  phases 
0.85  and  2.84,  i.e  phase  0.85  of  cycles  zero  and  two.  In  both  ex¬ 
amples,  the  /V-band  fluxes  are  lower  in  cycle  two  than  for  the 
same  phase  of  (a)  cycle  one  and  (b)  cycle  zero  by  up  to  ~100  Jy, 
corresponding  to  about  30%  or  one  to  two  standard  deviations. 
As  for  the  intracycle  photometry  variations  discussed  in  the  pre¬ 
vious  paragraph,  the  differences  are  most  pronounced  toward  the 
silicate  emission  feature  at  9.8  gm.  This  might  indicate  an  irreg¬ 
ularity  in  the  variability  cycle  in  the  A'-band.  In  the  optical  band 
the  visual  light  curve  clearly  shows  that  the  pulsation  cycles  are 
not  perfectly  symmetric  (see  Fig.  1). 

Summarizing,  our  data  exhibit  a  l-2<x  signature  of  intracy¬ 
cle,  as  well  as  cycle-to-cycle,  flux  variations  at  wavelengths  of 
8-1 3  gm,  which  are  most  pronounced  toward  the  silicate  emis¬ 
sion  feature  at  9.8y/m.  This  indication  is  consistent  with  obser¬ 
vations  by  Monnier  et  al.  (1998),  who  report  temporal  variations 
in  the  mid-infrared  spectra  of  late  type  stars.  In  particular,  stars 
with  a  strong  emission  feature  around  9.8  gm  showed  evident 
changes  in  the  spectral  profile.  Considerable  phase  variation  of 
the  /V-band  spectra  can  also  be  found  in  Alvarez  &  Plez  (1998) 
and  Guha  Niyogi  et  al.  (201 1). 


5.  Modeling  of  our  MIDI  data 

Self-consistent  models  describing  the  dynamic  atmosphere  of 
the  central  source  and  the  dust  shell  are  still  very  rare.  This  is 
particularly  true  for  oxygen-rich  stars.  Nevertheless,  some  ad¬ 
vance  in  this  domain  has  already  been  successfully  achieved 
(Ireland  &  Scholz  2006;  Hofner  &  Andersen  2007;  Ireland  et  al. 
2008;  Hofner  2008).  Here,  we  use  an  approach  of  an  ad  hoc 
radiative  transfer  modeling  of  the  dust  shell  where  the  central 
stellar  source  is  described  by  readily  available  and  established 
dust-free  dynamic  model  atmosphere  series,  as  introduced  by 
Wittkowski  et  al.  (2007). 

The  dust  shell  surrounding  the  central  star  is  modeled  by 
the  Monte  Carlo  radiative  transfer  code  mcsim_mpi  (Ohnaka 
et  al.  2006).  The  radiative  transfer  model  requires  an  assumption 
on  the  spectral  energy  distribution  (SED)  of  the  central  stellar 
source.  Since  we  expect  the  photosphere  and  molecular  layers 
of  RR  Aql  to  be  partly  resolved  with  our  MIDI  baselines  ( V 2  ~ 
0.4  for  0  =  10  mas,  B  =  120  m,  A  =  10/rm),  we  also  need  to 
model  the  intensity  distribution  across  the  atmospheric  layers. 


For  these  purposes,  we  use  the  dust-free  dynamic  model  atmo¬ 
spheres,  based  on  self-excited  pulsation  models  (P  and  M  se¬ 
ries),  by  Ireland  et  al.  (2004b, a,  and  references  therein)  as  the 
currently  best  available  option. 

The  P  and  M  series  were  constructed  to  reproduce  the 
M-type  Mira  prototypes  o  Cet  and  R  Leo.  The  hypothetical  non¬ 
pulsating  “parent”  stars  have  masses  M/M0  =  1 .0  (P  series),  and 
1.2  (M  series),  Rosseland  radius  Rp/RQ  =  241  (P)  and  260  (M), 
and  effective  temperatures  Teg  =  2860  K  (P)  and  2750  K  (M). 
The  models  have  solar  abundances,  luminosity  L/LQ  =  3470, 
and  a  pulsation  period  of  332  days,  close  to  the  periods  of  o  Cet 
(332  days)  and  R  Leo  (310  days).  They  are  available  for  20 
phases  in  three  cycles  for  the  M  series,  and  25  phases  in  four 
cycles  for  the  P  series.  These  models  have  been  shown  to  be 
consistent  with  broad-band  near-infrared  interferometric  data  of 
o  Cet  by  Woodruff  et  al.  (2004)  and  of  R  Leo  by  Fedele  et  al. 
(2005)  respectively.  They  have  also  been  shown  to  be  consis¬ 
tent  with  spectrally  resolved,  near-infrared  interferometric  data 
of  the  longer  period  (~430  days)  Mira  variable  S  Ori  exhibiting 
a  significant  variation  in  the  angular  size  as  a  function  of  wave¬ 
length,  which  was  interpreted  as  the  effect  of  molecular  layers 
lying  above  the  continuum-forming  layers  as  predicted  by  these 
model  series  (Wittkowski  et  al.  2008).  Tej  et  al.  (2003)  indicate 
that  the  spectral  shape  of  Mira  variables  can  be  reproduced  rea¬ 
sonably  well  even  for  stars  with  different  pulsation  periods.  In 
comparison  with  the  stellar  parameters  of  the  P  and  M  model  se¬ 
ries,  RR  Aql  pulsates  with  a  longer  period  of  394  days.  The  main 
sequence  precursor  mass  of  RR  Aql  is  1.2  ±  0.2  M/MQ  (Wyatt 
&  Cahn  1983),  and  its  spectral  type  is  M6e-M9  (Samus  et  al. 
2004),  versus  M5-M9  (P  series)  and  M6-M9.5  (M  series).  For 
the  purpose  of  the  modeling  effort  conducted  here,  we  can  as¬ 
sume  that  the  general  properties  of  the  P  and  M  model  series  are 
also  valid  for  a  longer  period  variable  such  as  RR  Aql.  However, 
an  exact  match  of  model  parameters  and  observed  parameters  as 
a  function  of  phase  may  not  be  expected. 

The  flux  and  visibility  values  of  the  combined  atmosphere 
and  dust  shell  (“global")  model  are  computed  for  each  spectral 
channel  as 


./total 

Ftotal 


/star  +  /dust 

Vstar  +  |^Fdust 
/total 


ratt 

./star 

./total 


(1) 

(2) 


where  f‘ ®  represents  the  attenuated  flux  from  the  dust-free 
model  atmosphere5  (extended  by  a  blackbody  approximation  for 
wavelengths  beyond  23  /urn),  /dust  represents  the  flux  from  the 
dust  shell  alone,  /totai  is  the  addition  of  these  two  terms,  and  Fstar 
and  Vdust  are  the  synthetic  visibilities  computed  for  the  dust-free 
model  atmosphere  and  for  the  dust  shell,  respectively. 


5. 1.  MIDI  model  parameters 

A  clear  silicate  feature  is  identified  in  the  spectra  of  RR  Aql 
(see  Fig.  3  and  Sect.  4).  RR  Aql  is  one  of  31  oxygen-rich  stars 
studied  by  Lorenz-Martins  &  Pompeia  (2000).  The  authors  mod¬ 
eled  the  stars  using  ALO3  and  silicate  grains  and  suggest  that 
the  dust  chemistry  of  RR  Aql  contains  only  silicate  grains.  In 
this  study  we  examined  both  of  these  two  dust  species,  ALO3 
grains  (Begemann  et  al.  1997;  Koike  et  al.  1995)  and  silicates 
grains  (Ossenkopf  et  al.  1992).  The  amount  of  dust  is  described 
by  the  optical  depths  at  To  =  0.55  pm.  The  grain  size  was  set 

5  Note  that  /“"  may  slightly  vary  across  the  stellar  disk,  which  is  not 
taken  into  account  here. 
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Table  2.  Best-fitting  model  parameters  for  each  epoch. 


Epoch 

®v,s 

Tv 

^in/  ^Phot 

P 

^Phot 

(sil.) 

(sil.) 

(sil.) 

[mas] 

A 

0.58 

3.0  ±  0.9 

4.5  ±  1.1 

2.5  ±  0.4 

6.6  ±  2.4 

B 

0.81 

3.0  ±  0.7 

4.5  ±  1.6 

2.5  ±  0.5 

7.1  ±  2.9 

C 

0.86 

2.5  ±  0.8 

5.0  ±  1.3 

2.5  ±  0.5 

8.0  ±2.8 

D 

1.53 

3.7  ±  0.6 

3.2  ±  1.1 

2.5  ±  0.5 

6.7  ±  2.6 

E 

1.77 

2.0  ±  0.8 

4.5  ±  1.4 

3.0  ±  0.4 

8.0  ±  2.9 

F 

2.46 

2.7  ±  0.8 

4.6  ±  1.1 

3.0  ±  0.5 

7.7  ±  2.4 

G 

2.55 

3.0  ±  0.6 

3.4  ±  1.7 

2.1  ±0.6 

7.2  ±  2.5 

H 

2.62 

2.0  ±  0.8 

4.5  ±  1.1 

2.4  ±  0.5 

8.5  ±  2.7 

I 

2.74 

2.0  ±  0.9 

4.0  ±  1.5 

2.6  ±  0.5 

8.5  ±  2.8 

J 

2.84 

3.2  ±  0.7 

3.0  ±  1.5 

2.0  ±  0.6 

7.2  ±  2.6 

K 

3.43 

5.0  ±0.8 

3.0  ±  1.3 

3.0  ±  0.5 

7.3  ±2.2 

L 

3.54 

2.3  ±  0.9 

5.0  ±  1.5 

2.5  ±  0.4 

8.2  ±  2.4 

M 

3.57 

2.5  ±  0.7 

4.5  ±  1.1 

2.5  ±  0.4 

7.5  ±  2.4 

Table  3.  Average  model  parameters. 

Model 

pv 

Tv 

^in/^Phot 

P 

^Phot 

(AI2O3) 

(silicate) 

(silicate) 

(silicate) 

[mas] 

M21n 

0.0 

2.8  ±0.8 

4.1  ±  0.7 

2.6  ±  0.3 

7.6  ±  0.6 

to  0.1  //m  for  all  grains.  We  set  the  photospheric  radius  to  the 
well-defined  continuum  photospheric  radius  at  A  =  1.04 /mi 
C^Phot  =  Wi.oi).  The  density  gradient  was  defined  by  a  single 
power  law  p(r)  o c  r~p  with  index  p.  The  shell  thickness  was  set 
tO  Rout/ R'm  =  1000. 

The  global  model  includes  seven  parameters.  The  radiative 
transfer  model  describing  the  circumstellar  dust  shell  includes 
six  parameters:  the  optical  depths  ry  (AI2O3)  and  ry(silicate), 
the  inner  boundary  radii  Rin/Rphot  (AI2O3)  and  Rin/Rphot  (sili¬ 
cate),  and  the  density  gradients  pA  (AI2O3)  and  pB  (silicate). 
The  model  of  the  P/M  atmosphere  series  is  another  parame¬ 
ter.  Here,  we  used  ten  models  covering  one  complete  cycle  of 
the  M  series:  M16n  (model  visual  phase  <Tmodei  =  0.60),  M18 
(0.75),  M18n  (0.84),  M19n  (0.90),  M20  (0.05),  M21n  (0.10), 
M22  (0.25),  M23n  (0.30),  M24n  (0.40),  and  M25n  (0.50). 

We  computed  a  grid  of  dust-shell  models  for  each  of  these  M 
models,  including  all  combinations  of  optical  depths  ry  (AI2O3) 
=  0.0,  0.1,  0.2,  0.5,  0.8;  tv  (silicate)  =  2.0,  2.5,  3.0,  3.5,  4.0, 

4.5,  5.0;  flin/flphot  (AI2O3)  =  2.0,  2.5,  3.0;  Rm/RPhot  (silicate)  = 

2.5,  3.5,  4.5,  5.5,  6.5;  pA  (A1203)  =  2.0,  2.5,  3.0,  3.5;  and  pB 
(silicate)  =  2.0,  2.5,  3.0,  3.5. 

In  a  first  selection  of  suitable  models,  we  compared  the  MIDI 
data  of  each  epoch  to  our  entire  grid  of  models.  The  angular  di¬ 
ameter  0phot  was  the  only  free  parameter.  We  increased  the  un¬ 
certainty  of  the  part  of  the  photometric  spectra  around  9.5  pm, 
which  is  strongly  affected  by  telluric  absorption.  Otherwise,  the 
weight  of  each  data  point  was  given  by  the  corresponding  uncer¬ 
tainty.  For  each  epoch  we  kept  the  best  -30  models  with  lowest 
X1  values.  In  a  next  step,  a  grid  with  finer  steps  around  these 
parameters  was  computed.  This  procedure  was  repeated  several 
times.  The  results  of  the  automatic  selection  were  visually  in¬ 
spected.  We  completed  the  selection  with  ten  models  for  each 
epoch  that  agrees  with  the  data  best. 

5.2.  MIDI  model  results 

For  each  of  the  13  epochs  we  found  the  best-fitting  model  pa¬ 
rameters  that  include  the  best-fitting  set  of  dust  parameters  of 
the  radiative  transfer  model,  together  with  the  best-fitting  model 
of  the  P/M  series.  We  used  the  procedure  outlined  in  Sect.  5.1. 


Following  Lorenz-Martins  &  Pompeia  (2000),  we  investi¬ 
gated  dust  shells  including  AFO3  and/or  silicate  shells  with  dif¬ 
ferent  inner  radii  and  density  gradients.  We  obtained  best-fit  re¬ 
sults  with  a  silicate  shell  alone,  and  the  addition  of  an  ALO3 
shell  did  not  result  in  any  improvement  in  the  model  fits.  This 
result  is  consistent  with  Lorenz-Martins  &  Pompeia  (2000),  who 
classify  RR  Aql  as  a  source  that  can  be  described  with  a  silicate 
dust  shell  alone. 

The  best-fitting  parameters  for  each  epoch  are  listed  in 
Table  2.  The  table  lists  the  epoch,  the  phase  at  the  epoch,  the  op¬ 
tical  depth  Ty,  the  inner  boundary  radius  Rin/Rphot,  the  density 
distribution  p,  and  the  continuum  photospheric  angular  diameter 
$Phot  -  Here,  the  dust  shell  parameters  are  those  of  the  silicate  dust 
shell.  The  errors  of  the  dust  shell  parameters  are  derived  in  the 
same  way  as  the  standard  deviation  based  on  the  find  best-fitting 
models  (Table  2).  Table  3  lists  the  average  model  parameters  of 
our  different  epochs.  The  average  phase  of  our  observations  is 
®y  =  0.64  ±  0.15. 

The  best-fitting  model  atmosphere  was  M2  In  (7nlo(|ci  = 
2550  K,  ®Modei=0.1)  f°r  all  epochs  covering  minimum  to  pre¬ 
maximum  pulsation  phases  (0.45-0.85).  The  difference  between 
the  average  phase  of  our  observation  and  the  phase  of  the  best¬ 
fitting  model  atmosphere  can  most  likely  be  explained  by  the 
different  stellar  parameters  of  RR  Aql  compared  to  those  of  the 
M  series.  It  means  that  Teg  of  RR  Aql  at  the  time  of  observation 
is  similar  to  Teg  of  model  M2  In,  which  has  a  different  phase,  but 
its  basic  parameters  are  also  not  exactly  those  of  RR  Aql. 

Figure  5  includes  the  model  flux  and  model  visibility  com¬ 
pared  to  the  observed  values  for  the  example  of  epoch  A.  For 
most  epochs,  the  agreement  between  the  models  and  the  ob¬ 
served  data  is  very  good,  in  particular  for  the  visibility  spectra. 
The  visibility  spectra  contain  information  about  the  radial  struc¬ 
ture  of  the  atmospheric  molecular  layers  and  the  surrounding 
dust  shells.  In  the  wavelength  range  from  8  pm  to  ~9pm,  the 
dust  is  fully  resolved  with  baselines  longer  than  -20  m,  and  the 
visibility  contribution  of  the  dust  reaches  into  the  second  lobe 
of  the  visibility  function.  RR  Aql  is  characterized  by  a  partially 
resolved  stellar  disk  that  includes  atmospheric  layers,  with  a  typ¬ 
ical  drop  in  the  visibility  function  -  10//m,  where  the  flux  contri¬ 
bution  of  the  silicate  emission  is  highest  and  the  flux  contribution 
of  the  star  relative  to  the  total  flux  decreases.  Beyond  ~10//m, 
the  spatially  resolved  radiation  from  the  optically  thin  dust  shell 
starts  to  be  a  notable  part  of  the  observed  total  flux.  From  -10  pm 
to  13//m,  the  dust  contribution  becomes  nearly  constant  while 
the  stellar  contribution  increases,  and  this  results  in  a  rebound  of 
the  visibility  function.  Photometry  spectra  also  fit  well.  For  some 
epochs,  the  models  predict  higher  fluxes  than  observed  near  the 
silicate  emission  feature  at  -9.8  pm.  Generally,  our  attempt  at  a 
radiative  transfer  model  of  the  circumstellar  dust  shell  that  uses 
dynamical  model  atmosphere  series  to  describe  the  central  stel¬ 
lar  source  can  reproduce  the  shape  of  both  the  visibility  and  the 
photometry  spectra  very  well.  Our  model  also  follows  the  shape 
of  the  SED  in  the  range  of  l-40//m  (Fig.  3). 

The  obtained  model  parameters  for  the  different  epochs  do 
not  indicate  any  significant  dependence  on  phase  or  cycle.  This  is 
consistent  with  the  result  from  Sect.  4. 1  that  a  direct  comparison 
of  visibility  values  of  different  phases  or  cycles  did  not  show 
any  variability  within  our  uncertainties  and  that  the  photometric 
values  indicated  only  small  variations  of  up  to  -2 cr.  However, 
we  cannot  confirm  or  deny  an  eventual  phase  dependence  of  the 
dust  formation  process  that  affects  the  8-1 3  pm  visibility  and 
photometry  values  by  less  than  our  uncertainties  of  about  5-20% 
and  10-50%,  respectively. 
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A  silicate  dust  shell  alone,  i.e.  without  the  addition  of  an 
AI2O3  dust  shell,  provides  the  best  agreement  with  our  data.  The 
average  optical  depth  of  the  silicate  dust  is  ry(silicate)  =  2. 8+0. 8 
at  A  =  0.55  pm  (corresponding  to  0.03  at  A  =  8  pm.  0.06  at 
A  =  12;um,  and  a  maximum  within  8-12/rm  of  0.22  at  A  - 
9.8pm).  The  inner  radius  of  the  dust  shell  expressed  in  Rphot  at 
A  —  1.04pm  is  Rin  =4.1+  0.7/?phot.  The  intensity  profile  of 
the  dust-free  dynamic  model  atmosphere  extends  to  about  1.5- 
2  photosphere  radii  at  10pm  (cf.  Wittkowski  et  al.  2007).  The 
power-law  index  of  the  density  distribution  is  p  =  2. 6+0. 3.  Wind 
models  predict  a  power-law  index  of  the  density  distribution  of 
2  at  radii  outward  the  dust  formation  zone.  However,  for  radii 
close  to  the  dust  formation  zone  ( r  <~  10  7?*),  a  larger  index 
is  expected  (cf,  e.g.,  Fig.  6  in  Wittkowski  et  al.  2007),  which  is 
consistent  with  our  result. 

The  average  photospheric  angular  diameter  results  in  6>phot  = 
7.6  ±  0.6  mas.  This  value  is  lower  than  the  A'-band  (A  =  2.2pm, 
AT  =  0.4pm)  UD  diameter  of  $ud  =  10.73  ±  0.66  mas  de¬ 
rived  by  van  Belle  et  al.  (2002)  at  a  minimum  phase  of  0.48. 
This  can  most  likely  be  explained  by  the  different  radius  def¬ 
initions.  Our  diameter  is  a  photospheric  angular  diameter  that 
has  already  been  corrected  for  effects  of  molecular  layers  ly¬ 
ing  above  the  continuum  photosphere  using  the  prediction  by 
the  best-fitting  model  atmosphere,  while  the  diameter  by  van 
Belle  et  al.  is  a  broad-band  uniform  disk  diameter  that  still  con¬ 
tains  contamination  by  molecular  layers.  It  is  known  that  broad¬ 
band  UD  diameters  may  widely  overestimate  the  photospheric 
radius  owing  to  contamination  by  molecular  layers  lying  above 
the  photosphere  (cf.,  e.g.,  Ireland  et  al.  2004b;  Fedele  et  al.  2005; 
Wittkowski  et  al.  2007).  With  the  parallax  of  n  =  1.58  +  0.40  mas 
(Vlemmings  &  van  Langevelde  2007),  our  value  for  0phot  corre¬ 
sponds  to  a  photospheric  radius  of  Rpgot  =  520  '  R0.  Together 

with  the  bolometric  magnitude  of  RR  Aql  of  //fhoi  =  3.71 
and  Amboi  =  1.17  from  Whitelock  et  al.  (2000),  our  value  for 
Ophot  at  <b v  ~  0.64  corresponds  to  an  effective  temperature  of 
Tefi  ~  2420  ±  200  K.  This  value  is  consistent  with  the  effective 
temperature  of  the  best-fitting  model  atmosphere  M2  In,  which 
is  Teg  =  2550  K. 


5.3.  Model  simulations 

We  performed  a  number  of  model  simulations  in  order  to  inves¬ 
tigate  the  visibility  and  photometry  variations  that  are  theoret¬ 
ically  expected  in  the  8- 13  pm  wavelength  range  for  the  typi¬ 
cal  parameters  of  RR  Aql  as  determined  above  in  Sect.  5.2.  We 
based  our  simulations  on  a  global  model  with  typical  parameters 
similar  to  those  derived  for  RR  Aql  (cf.  Table  3),  and  investigated 
the  effects  of  expected  variations  in  the  atmosphere  model  and 
the  dust  shell  parameters  during  a  pulsation  cycle  on  the  observ¬ 
able  photometry  and  visibility  spectra.  We  mainly  used  a  mean 
projected  baseline  length  of  60  m.  As  in  Sect.  5,  the  global  model 
consists  of  a  radiative  transfer  models  of  the  AEO3  and  silicate 
dust  shells  where  the  central  source  is  described  by  a  dust-free 
dynamic  model  atmosphere.  Table  4  lists  the  performed  simu¬ 
lations  including  the  phase  of  the  model,  the  optical  depth 
(AI2O3)  and  Ty  (silicate),  the  inner  boundary  radii  of  the  dust 
shell  Rm/ Rphot  (AI2O3)  and  R;n/Rphot  (silicate),  the  power-law 
indices  of  the  density  distributions  p  (AFO?)  and  p  (silicate), 
and  the  continuum  photospheric  angular  diameter  9pilol.  We  as¬ 
sume  a  phase  dependence  of  the  angular  photospheric  diameter 
of  about  20%  (Thompson  et  al.  2002a;  Ireland  et  al.  2004b). 
Figures  10-14  show  the  results  from  the  simulations.  Each 
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Fig.  10.  Simulation  1.  Synthetic  flux  (top)  and  visibility  amplitude  ( bot¬ 
tom )  in  the  wavelength  range  8-13  pm.  The  solid  lines  represent  the 
global  models,  the  dashed-dotted  lines  denote  the  unattenuated  stellar 
contribution  i.e.  Vstar  in  Eq.  (2),  the  dotted  lines  denote  the  attenuated 
stellar  contribution,  and  the  dashed  lines  denote  the  dust  shell  contribu¬ 
tions.  This  simulation  compares  two  models  consisting  of  the  same  dust 
shell  parameters  but  different  model  atmospheres,  the  post-maximum 
model  atmosphere  M21n  (4>m0dei  =  0.10),  and  the  minimum  model  at¬ 
mosphere  M23n  (4>m0dei  =  0.30).  The  photospheric  angular  diameter  is 
assumed  to  be  larger  at  post-maximum  phase  than  at  minimum  phase 
(8.5  mas  compared  to  7  mas).  The  projected  baseline  length  is  60  m. 
For  the  exact  model  parameters,  see  Table  4.  They  describe  variations 
around  our  best-fitting  model  for  RR  Aql. 

simulation  compared  two  models,  where  one  or  more  of  the 
model  parameters  have  been  varied. 

Simulation  1  (Fig.  10)  In  simulation  1,  we  investigate  the  effect 
of  model  atmospheres  at  different  pulsation  phases  and  effective 
temperatures.  We  compared  model  M21n  at  post-maximum  pul¬ 
sation  phase  (Omodei  =  0.10,  Teg  -  2550  K)  and  model  M23n  at 
pre-minimum  pulsation  phase  (<hmodei  =  0.30,  7’cii  =  2230  K), 
applying  the  same  dust  shell  parameters  for  both  models.  Here, 
the  phase  difference  between  M21n  and  M23n  of  0.2  corre¬ 
sponds  to  the  phase  difference  of  our  observations  (-0.64  ±  0.2) 
The  dust  parameters  were  based  on  the  average  parameters  (see 
Tab.  3)  of  RR  Aql  derived  in  Sect.  5.2.  The  projected  baseline 
length  was  set  to  Bp  =  60  m,  which  is  a  typical  mean  value 
of  our  MIDI  observations.  The  overall  shapes  of  the  visibil¬ 
ity  spectra  Vtotai  (M2 In)  and  Vtotai  (M23n)  are  similar.  The  dif¬ 
ferences  between  the  models  are  wavelength-dependent  in  the 
range  of  5-20%.  The  shape  of  the  photometry  spectra  are  sim¬ 
ilar  for  both  tested  models  /tota  1  (M21n)  and  /totai  (M23n).  The 
largest  difference  of  -25%  is  seen  around  1 0 pm,  with  less  flux 
at  minimum  pulsation  phase.  Compared  to  the  uncertainties  of 
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Table  4.  Model  simulations,  comparing  two  models  that  differ  in  one  or  more  parameters  (marked  by  bold  face). 


Simulation 

Model 

®vis 

Tv 

A1(203) 

tva 

Tv 

(silicate) 

tvb 

^in/  ^Phot 

(A1203) 

ra 

Rin  /  Rpiiol 

(silicate) 

rb 

P 

(A1203) 

pa 

P 

(silicate) 

pb 

llpilul 

[mas] 

diam. 

Projected 
baseline  [m] 
Bp 

1 

M21n 

0.1 

0.0 

3.0 

2.2 

4.5 

2.5 

2.5 

8.5 

60 

M23n 

0.3 

0.0 

3.0 

2.2 

4.5 

2.5 

2.5 

7.0 

60 

2 

M21n 

0.1 

0.0 

3.0 

2.2 

4.5 

2.5 

2.5 

8.5 

60 

M23n 

0.3 

0.0 

4.0 

2.2 

3.8 

2.5 

2.5 

7.0 

60 

3a 

M21n 

0.1 

0.0 

3.0 

2.2 

4.5 

2.5 

2.5 

8.5 

60 

M21n 

0.1 

0.0 

5.0 

2.2 

3.0 

2.5 

2.5 

8.5 

60 

3b 

M21n 

0.1 

0.0 

3.0 

2.2 

4.5 

2.5 

2.5 

8.5 

40 

M21n 

0.1 

0.0 

5.0 

2.2 

3.0 

2.5 

2.5 

8.5 

40 

3c 

M21n 

0.1 

0.0 

3.0 

2.2 

4.5 

2.5 

2.5 

8.5 

25 

M21n 

0.1 

0.0 

5.0 

2.2 

3.0 

2.5 

2.5 

8.5 

25 

3d 

M21n 

0.1 

0.0 

3.0 

2.2 

4.5 

2.5 

2.5 

8.5 

15 

M21n 

0.1 

0.0 

5.0 

2.2 

3.0 

2.5 

2.5 

8.5 

15 

4 

M21n 

0.1 

0.0 

3.0 

2.2 

4.5 

2.5 

2.5 

8.5 

60 

M21n 

0.1 

0.6 

3.0 

2.2 

4.5 

2.5 

2.5 

8.5 

60 

Notes.  Here,  ®vis  is  the  phase,  tv  the  optical  depths,  Rin/Rphot  the  inner  boundary  radii,  p  the  density  gradients,  0Phot  the  angular  diameter. 


Fig.  11.  Simulation  2.  Like  simulation  1  in  Fig.  10,  but  for  the  parame¬ 
ters  of  simulation  2  comparing  a  post-maximum  atmosphere  model  and 
a  minimum  model  as  in  simulation  1,  but  where  the  dust  is  also  assumed 
to  be  closer  to  the  star  with  larger  optical  depth  at  minimum  phase  and 
farther  from  the  star  with  lower  optical  depth  at  post-maximum  phase. 


our  MIDI  measurements  of  ~5-20%  in  the  visibility  spectra  and 
of  ~  20%  on  average  in  the  photometric  spectra,  simulation  1  is 
consistent  with  our  non-detection  of  visibility  variations  and  a 
marginal  detection  of  photometry  variations. 


Fig.  12.  Simulation  3.  As  for  simulation  2  in  Fig.  11,  but  where  the 
atmosphere  model  is  not  varied. 


Simulation  2  (Fig.  11)  Compared  to  simulation  1,  we  adjusted 
the  parameters  according  to  the  assumption  that  the  outer  lay¬ 
ers  of  the  atmosphere  are  cooler  near  minimum  visual  pulsa¬ 
tion  phase,  and  therefore  more  dust  grains  can  form,  and  a 
higher  molecular  opacity  can  be  expected.  This  theoretical  as¬ 
sumption  is  consistent  with  mid-infrared  interferometric  obser¬ 
vations  (Wittkowski  et  al.  2007),  where  the  observed  data  indi¬ 
cated  more  dust  formation  near  minimum  pulsation  phase,  with 
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Wavelength  [/xm] 

Fig.  13.  Simulation  3.  As  Fig.  12,  but  showing  the  visibility  spectra  for 
different  projected  baseline  lengths  of  (from  top  to  bottom)  40  m,  25  m, 
and  15  m. 


inner  boundary  radii  of  the  dust  shell  located  closer  to  the  star. 
We  increased  the  optical  depth  for  model  M23n.  We  also  set  a 
more  extended  dust  shell  after  visual  maximum  for  model  M2  In 
in  comparison  to  the  model  M23n.  This  setting  results  in  spec¬ 
tra  with  differences  in  the  range  of  2-19%  for  the  visibility  val¬ 
ues.  The  photometry  differs  again  mostly  around  10  pm  with  a 
maximum  difference  of  20%.  Compared  to  simulation  1,  the  re¬ 
sults  are  very  similar,  but  with  less  differences  in  the  photometry 
spectra.  As  for  simulation  1,  the  results  are  consistent  with  our 
not  detecting  visibility  variations  and  marginally  detecting  pho¬ 
tometry  variations  within  our  measurement  uncertainties. 

Simulation  3  (Figs.  12  and  13)  Here,  we  investigate  the  effect  of 
different  dust  shell  parameters  (as  in  simulation  2),  but  keeping 
the  M  model  constant.  We  used  model  M21n  (phase  0.1),  which 
was  the  best-fitting  model  to  our  RR  Aql  epochs  and  also  used  a 
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Fig.  14.  Simulation  4.  Comparison  of  a  model  with  a  silicate  shell  only 
and  a  model  that  also  includes  an  AI1O3  dust  shell  with  lower  optical 
depth.  For  the  model  parameters,  see  Table  4;  For  the  description  of  the 
figure,  see  the  caption  of  Fig.  10. 

constant  photospheric  angular  diameter.  We  compared  two  (sili¬ 
cate)  dust  models  as  in  simulation  2,  where  in  one  model  the  op¬ 
tical  depth  is  lower  and  the  inner  radius  of  the  dust  shell  larger 
compared  to  the  other  model.  These  sets  of  model  parameters 
lead  to  very  similar  results,  where  the  visibility  spectra  almost 
overlap,  with  a  maximum  difference  of  4%  for  a  projected  base¬ 
line  length  of  60  m  (Fig.  12).  The  photometry  spectra  are  also 
very  similar  with  a  maximum  difference  of  7%.  Figure  1 3  shows 
the  visibility  results  for  different  projected  baseline  lengths  of 
40  m,  25  m,  and  15  m.  With  projected  baseline  lengths  of  less 
than  60  m,  the  difference  between  the  visibility  values  become 
more  pronounced  (10%  difference  beyond  1 0 pm  with  projected 
baseline  lengths  of  40  m  and  25  m).  This  can  be  understood  be¬ 
cause  the  dust  shell  is  over-resolved  with  a  projected  baseline 
length  of  60  m,  so  that  the  visibility  cannot  be  sensitive  to  vari¬ 
ations  in  the  geometry  of  the  dust  shell.  With  lower  projected 
baseline  length,  the  visibility  values  are  located  in  the  first  lobe 
of  the  visibility  function,  hence  more  sensitive  to  variations  of 
the  dust  shell  geometry.  However,  with  the  even  lower  projected 
baseline  lengths  of  15  m,  the  difference  between  individual  mod¬ 
els  again  starts  to  be  less  pronounced  (<7%).  This  indicates  an 
optimum  projected  baseline  length  of  ~25-40  m  to  characterize 
the  dust  shell  geometry  for  sizes  as  found  for  RR  Aql. 

Simulation  4  (Fig.  14)  Furthermore,  we  investigate  the  effect  of 
adding  small  amounts  of  AFO3  grains  to  a  thus  far  pure  silicate 
dust  shell.  We  examined  ryfAhO?)  of  0.3  and  0.6.  Figure  14 
shows  the  simulation  with  rytAFCF)  of  0.6,  showing  that  the 
addition  of  small  amounts  of  AFO?  grains  leads  to  almost  iden¬ 
tical  photometry  and  visibility  spectra  (differences  <2%). 
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6.  Discussion 

Our  MIDI  observations  of  RR  Aql  did  not  show  significant  vari¬ 
ations  in  the  8-1 3  /urn  visibility  within  the  examined  pulsation 
phases  between  0.45  and  0.85  within  a  total  of  four  pulsation  cy¬ 
cles,  and  showed  only  marginal  variations  in  the  8-1 3  /urn  flux. 

We  performed  model  simulations  with  expected  variations 
of  the  pulsation  phase  of  the  innermost  dust-free  atmosphere 
model  and  of  the  parameters  of  the  surrounding  dust  shells. 
These  model  simulations  show  that  visibility  variations  are  in¬ 
deed  not  expected  for  the  parameters  and  observational  settings 
of  RR  Aql  at  wavelengths  of  8-13/rm  within  the  uncertainties 
of  our  observations.  Variations  in  the  flux  spectra  may  in  some 
cases  just  be  detectable.  Thus,  our  observational  result  of  a  con¬ 
stant  visibility  and  only  slightly  varying  flux  at  wavelengths  of 
8-13  /jm  are  consistent  with,  and  not  contradicting,  theoretical 
expectations  of  a  pulsating  atmosphere. 

Our  model  simulations  indicate  that  detections  of  pulsation 
effects  at  mid-infrared  wavelengths  would  in  particular  benefit 
from  smaller  uncertainties  in  the  photometric  spectrum  than  in 
our  current  data.  Also,  a  wide  range  of  projected  baselines  at 
each  phase,  for  RR  Aql  in  particular  including  baseline  lengths 
around  20-30  m,  would  help  us  to  distinguish  models  of  different 
pulsation  phases. 

For  our  analysis,  we  have  considered  silicate  (Ossenkopf 
et  al.  1992)  and  AhO?  (Begemann  et  al.  1997;  Koike  et  al. 
1995)  dust  species,  following  the  work  by  Lorenz-Martins  & 
Pompeia  (2000),  and  used  fixed  grain  sizes  of  0. 1  /um.  It  is  known 
that  other  dust  species  with  a  more  complex  gran  size  distribu¬ 
tions  may  occur  in  the  circumstellar  environment  of  Mira  vari¬ 
ables  (e.g.,  Hofner  2008;  Molster  &  Waters  2003,  and  references 
therein).  We  showed  that  a  model  including  only  a  silicate  dust 
shell  can  reproduce  the  observed  RR  Aql  visibility  and  flux  spec¬ 
trum  at  8-1 3 /-I m.  The  addition  of  an  AFO3  dust  shell  with  com¬ 
parable  low  optical  depth  did  not  significantly  improve  the  fit  to 
our  data.  However,  our  model  simulations  have  shown  that  our 
8-13  fjm  visibility  and  flux  values  are  not  sensitive  to  the  addi¬ 
tion  of  an  AI1O3  dust  shell  with  low  optical  depth  within  our  un¬ 
certainties.  As  a  result,  we  cannot  exclude  the  presence  of  an  in¬ 
ner  AI2O3  dust  shell  in  addition  to  the  silicate  dust  shell.  Woitke 
(2006)  uses  dynamical  models  for  dust-driven  winds  of  oxygen- 
rich  AGB  stars,  including  frequency-dependent  radiative  trans¬ 
fer,  and  finds  that  dust  temperatures  strongly  depend  on  material. 
Two  dust  layers  are  formed  in  his  dynamical  models,  almost  pure 
glassy  AhO?  close  to  the  star  (r  >1.5  Rs tar)  and  the  more  opaque 
Fe-poor  Mg-Fe-silicates  farther  out  at  4-5  Rsta r. 

The  dust-free  dynamic  model  atmospheres  predict  a  signif¬ 
icant  dependence  of  the  characteristics  of  the  molecular  lay¬ 
ers  on  the  stellar  pulsation  phase  at  near-infrared  wavelengths 
(Ireland  et  al.  2004b, a).  There  are  a  few  near-infrared  interfer¬ 
ometric  observations,  which  detected  a  clear  variation  in  the 
continuum  angular  diameters  with  pulsation  phase  (Thompson 
et  al.  2002b;  Perrin  et  al.  1999;  Young  et  al.  2000b;  Woodruff 
et  al.  2004).  Spectrally  resolved  near-infrared  interferometric 
measurements  at  different  phases,  such  as  the  AMBER  obser¬ 
vations  by  Wittkowski  et  al.  (2008)  but  at  more  than  one  phase, 
promise  to  lead  to  stronger  constraints  on  dynamic  model  atmo¬ 
spheres. 

Our  model  simulations  based  on  the  combination  of  dy¬ 
namic  model  atmospheres  with  a  radiative  transfer  model  of 
the  dust  shell  predict  only  small  variations  with  phase  at  mid- 
infrared  wavelengths.  The  difference  near  a  wavelength  of  10  gm 
amounts  to  ~5%  for  the  visibility  values  and  to  ~25%  for  the 
photometry  values.  These  results  lead  to  the  suggestion  that 


the  stellar  photosphere  and  overlying  molecular  layers  pulsate, 
which  is  demonstrated  by  the  diameter  variations  in  the  near- 
infrared,  but  that  these  pulsations  cannot  be  detected  for  RR  Aql 
by  our  8- 13 /mi  interferometry  within  our  uncertainties.  In  addi¬ 
tion,  our  visibility  uncertainties  do  not  allow  us  to  exclude  vari¬ 
ations  in  the  geometry  and  optical  depth  of  the  dust  shell  as  a 
function  of  pulsation  phase  as  observed  previously  for  other  tar¬ 
gets  (Lopez  et  al.  1997).  Possible  explanations  for  not  detecting  a 
phase-dependence  of  the  dust  shell  parameters  in  our  study  com¬ 
pared  to  the  ISI  observations  include:  (i)  longer  baseline  lengths 
in  our  study  that  often  fully  resolve  the  very  extended  silicate 
dust  shell;  (ii)  our  limited  phase  coverage  between  minimum  and 
pre-maximum  (0.45-0.85)  phases.  Possible  variations  over  the 
whole  pulsation  cycle  cannot  be  excluded. 

The  observed  variability  of  the  8-1 3  /mi  flux  at  a  significance 
level  of  1— 2cr  may  indicate  variations  in  the  stellar  radiation  re¬ 
emitted  by  the  dust  and/or  changes  in  either  the  dust  geometry 
or  optical  depth. 


7.  Summary  and  conclusions 

We  have  investigated  the  circumstellar  dust  shell  and  character¬ 
istics  of  the  atmosphere  of  the  oxygen-rich  Mira  variable  RR  Aql 
using  mid-infrared  interferometric  observations.  We  observed 
RR  Aql  with  the  VLTI/MIDI  instrument  at  different  pulsation 
phases  in  order  to  monitor  the  photometry  and  visibility  spectra. 
A  total  of  57  observations  were  combined  into  13  epochs  cover¬ 
ing  four  pulsation  cycles  between  April  2004  and  July  2007,  and 
covering  pulsation  phases  between  minimum  and  pre-maximum 
phases  (0.45-0.85). 

We  modeled  the  observed  data  with  an  ad-hoc  radiative 
transfer  model  of  the  dust  shell  using  the  radiative  transfer 
code  mcsim_mpi  by  Ohnaka  et  al.  (2006).  In  this  way,  we  used 
a  series  of  dust-free  dynamic  model  atmospheres  based  on 
self-excited  pulsation  models  (M  series,  Ireland  et  al.  2004b, a) 
to  describe  the  intensity  profile  of  the  central  source.  This 
study  represents  the  first  comparison  between  interferometric 
observations  and  theoretical  models  over  an  extended  range  of 
pulsation  phases  covering  several  cycles. 

Our  main  observational  results  are  as  follows: 

-  the  interferometric  data  do  not  show  any  evidence  of  intra¬ 
cycle  visibility  variations; 

-  the  data  do  not  show  any  evidence  of  cycle-to-cycle  visibility 
variations; 

-  The  8-13  /jm  flux  suggests  intracycle  and  cycle-to-cycle 
photometry  variations  at  a  significance  level  of  l-2cr. 
Follow-up  observations  with  higher  accuracy  using  a  dedi¬ 
cated  photometric  instrument,  such  as  VISIR  at  the  VLT,  are 
needed  to  confirm  this  result. 

These  observational  results  can  be  explained  by  dynamic 
model  atmospheres  and  variations  in  the  dust  shell  parameters. 
Simulations  using  different  phases  of  the  dynamic  model  atmo¬ 
sphere  and  different  sets  of  dust  shell  parameters  predict  vis¬ 
ibility  variations  that  are  lower  than  or  close  to  the  observed 
visibility  uncertainties  (5-20%).  Model-predicted  variations  of 
the  photometry  spectra  are  largest  around  wavelengths  of  10  /mi 
with  difference  of  up  to  ~25%  corresponding  to  up  to  l-2cr. 

The  best-fitting  model  for  our  average  pulsation  phase  of 
<T> v  =  0.64  ±  0.15  includes  a  silicate  dust  shell  with  an  op¬ 
tical  depth  of  Tv(silicate)  =  2.8  ±  0.8,  an  inner  radius  of 
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Rin  —  4.1  ±  0.7  flphot,  and  a  power-law  index  of  the  density  dis¬ 
tribution  of  p  =  2.6  ±  0.3.  The  corresponding  best-fitting  at¬ 
mosphere  model  of  the  series  used  to  describe  the  central  in¬ 
tensity  profile  is  M21n  (rm0dei  =  2550K,  d>m0dei  =  0.1)  with 
a  photospheric  angular  diameter  of  0phot  =  7.6  ±  0.6  mas.  The 
photospheric  angular  diameter  corresponds  to  a  photospheric 
radius  of  Rp hot  =  520^2  R0  and  an  effective  temperature  of 
Teg  ~  2420  ±  200  K.  The  latter  value  is  consistent  with  the 
effective  temperature  of  the  used  model  M2  In.  The  combined 
model  can  reproduce  the  shape  and  features  of  the  observed  pho¬ 
tometry  and  visibility  spectra  of  RR  Aql  well,  as  well  as  the  SED 
at  l-40jum. 

We  conclude  that  our  RR  Aql  data  can  be  described  by  a  sil¬ 
icate  dust  shell  surrounding  a  pulsating  atmosphere,  consistent 
with  observations  of  Mira  variables  (Lorenz-Martins  &  Pompeia 
2000).  The  effects  of  the  pulsation  on  the  mid-infrared  flux  and 
visibility  values  are  expected  to  be  less  than  about  25%  and  20%, 
respectively,  and  are  too  low  to  be  detected  within  our  measure¬ 
ment  uncertainties.  Although  the  addition  of  an  AEO3  dust  shell 
did  not  improve  the  model  fit,  our  simulations  also  indicate  that 
we  cannot  exclude  the  presence  of  an  inner  APO  ?  dust  shell  with 
relatively  low  optical  depth,  which  may  be  an  important  contrib¬ 
utor  to  the  dust  condensation  sequence. 

Our  modeling  attempt  at  a  radiative  transfer  model  of  the 
dust  shell  surrounding  a  dynamic  dust-free  model  atmosphere 
provides  constraints  on  the  geometric  extensions  of  atmospheric 
molecular  and  dust  shells.  It  should  be  noted  that  it  cannot  ex¬ 
plain  the  mechanism  by  which  the  observed  mass  loss  is  pro¬ 
duced  and  the  wind  is  driven.  Newer  models  of  M  type,  i.e.  of 
oxygen-rich  Miras,  that  include  atmospheric  dust  (e.g.,  Ireland 
et  al.  2008,  Ireland,  in  prep.)  typically  predict  effective  accel¬ 
eration  very  close  to  zero  in  high  layers  but,  so  far,  no  clear 
outward  acceleration,  i.e.  no  wind.  Subtle  details  are  being  dis¬ 
cussed  that  may  overcome  this  shortcoming  of  dynamic  mod¬ 
els  of  M  type  Miras  (e.g.,  Hofner  2011).  Future  observations 
aiming  at  characterizing  and  constraining  other  new  models  of 
the  mass-loss  process  and  the  wind  driving  mechanism  at  mid- 
infrared  wavelengths  would  benefit  from  obtaining  more  precise 
photometry  values  using  a  dedicated  instrument  like  VISIR  at 
the  VLT,  the  addition  of  shorter  baselines  to  characterize  the  ex¬ 
tension  of  the  silicate  dust  shell,  and  a  more  complete  coverage 
of  the  pulsation  cycle.  The  addition  of  concurrent,  spectrally  re¬ 
solved,  near-infrared  interferometry  would  be  needed  to  more 
strongly  constrain  atmospheric  molecular  layers  located  close  to 
the  photosphere. 
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Fig.  5.  VLTI/MIDI  interferometry  at  8-13  fim  of  RR  Aql  for  the  example  of  epoch  A  (stellar  phase  0.58)  and  epoch  I  (stellar  phase  2.75)  .  The 
panels  show  a)  the  flux;  b)  the  visibility  amplitude;  c)  the  corresponding  UD  diameter;  d)  the  corresponding  Gaussian  FWHM  diameter  as  a 
function  of  wavelength.  The  gray  shading  indicates  the  wavelength  region  around  9.5  /mi  that  is  affected  by  atmospheric  absorption.  Panel  e) 
shows  the  visibility  amplitude  as  a  function  of  spatial  frequency  for  three  averaged  bandpasses  of  8-9  /mi,  10-11  /mi,  and  12-13  /mi.  The  crosses 
with  error  bars  denote  the  measured  values.  The  solid  lines  indicate  our  best-fitting  model,  as  described  in  Sect.  5.  It  consists  of  a  combination 
of  a  dust-free  dynamic  model  atmosphere  representing  the  central  star  and  a  radiative  transfer  model  representing  the  surrounding  dust  shell.  The 
contributions  of  the  stellar  and  dust  components  alone  are  indicated  by  the  dotted  and  the  dashed  lines,  respectively. 
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